Decoherence versus disentanglement for two qubits in a squeezed bath. 
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We study the relation between the sudden death and revival of the entanglement of two qubits 
in a common squeezed reservoir, and the normal decoherence, by getting closer to the Decoherence 
Free Subspace and calculating the effect on the death and revival times. 
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I. INTRODUCTION 



In one-party quantum systems, coherence is destroyed 
by the action of the environment, a phenomena that is 
local and occurs asymptotically in time. 

On the other hand, there could be multiparty systems, 
with non-local quantum correlations, often referred to as 
quantum entanglement. The non-locality and coherence 
of the quantum entangled states makes them very im- 
portant in applications such as quantum teleportation 

[I] , quantum cryptography [2], dense coding [3] , etc. 

If the environment would act on the various parties the 
same way it acts on single systems, one would expect that 
a measure of entanglement, say the concurrence, would 
also decay exponentially in time. 

However, this is not always the case. Recently Yu and 
Eberly [4-6] showed that under certain conditions, the 
dynamics could be completely different and the quantum 
entanglement may vanish in a finite time. They called 
this effect "entanglement sudden death". This effect has 
also been observed experimentally [7]. 

A popular measure of entanglement is the concurrence 

[II] . For pure states is defined as: 

C(vl/) = |(*|$)|, (1) 

where: 

|^)i2 = (cr^ (8) (jI) (8) |^)i2, for two qubits. (2) 

For a mixed state, the concurrence is defined as: 

C{p) = max{0, y/X~i-y^-y/X^- VaI}, (3) 

where the are the eigenvalues (Ai being the largest 
one) of a non-Hermit ian matrix pp, and p is defined as: 

P={(Tl^al)p%al0al), (4) 

p* being the complex conjugate of p. 

In this paper we will explore the relation between the 
sudden death (and revival) of the entanglement between 
the two two-level atoms in a squeezed bath and the 
normal decoherence and the decoherence free subspace 
(DFS), which in this case is a two-dimensional plane. 



II. THE MODEL 

In the present work, we consider two two-level atoms 
that interact with a common squeezed reservoir, and we 
will focus on the evolution of the entanglement between 
them, using as a basis, the Decoherence Free Subspace 
states, as defined in references [8] and [9]. 

The master equation, in the Interaction Picture, for a 
two- level system in a broadband squeezed vacuum bath 
is given by [10]: 

d 1 

^ = + l)(2crpcr'^ - crVp - pcrV) 

-h ^'yN{2a^ pa — aa^ p — paa^) 

- ^7Me^^(2crV^T^ - crVV - pcrVt) 

— ^7Me*^ (2(7^(7 — aap — paa), (5) 

where 7 is th e spontane ous emission rate and N = 
sinh^ r^M = y^^]V(7VH-T) and ^ are the squeeze param- 
eters of the bath and a''', aare the usual Pauli raising and 
lowering matrices. 

It is simple to show that the above master equation 
can also be written in the Lindblad form with a single 
Lindblad operator S 

^ = ^^{2SpS^-S^Sp-pS^S), (6) 

with 

S = V]vTT(a) - ^e^^(at). 

For a two two-level system, the master equation has 
the same structure, but now the S operator becomes: 

S = VNTl{ai + ^2) - VNe'^{al + al) 

= cosh(r)(cri -h (J2) — sinh(r)e*^(cri + (tI)- (7) 

The Decoherence Free Subspace was found in ref [8] 
and consists of the eigenstates of S with zero eigenvalue. 
The states defined in this way, form a two-dimensional 
plane in Hilbert Space. Two orthogonal vectors in this 
plane are: 
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We can also define the states l^s) and |04) orthogonal 
to the {|0i) , 102)} plane: 



l<^4) 
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ViV2 + M2 



(M| + +)-7Ve-^^|--)). (11) 



To solve the master equation, we are going to use the 
basis {|0i), 102), l^^s)^ 1^4)} This solution depends on the 
initial state. We present the general solution in the Ap- 
pendix. 

In general, for density matrices written in the standard 
basis of the form: 



Pit) = 



/pii pi4\ 

P22 P23 

P32 P33 

\P41 P44/ 



(12) 



one easily finds [12, 13] that the concurrence is given 
by: C{p) = max{0,Cl{p),C2{p)}, where: 




FIG. 1: Concurrence as function of N, for \(pi) as the initial 
state. 



b) If we now consider 102) as the initial state, we get 
the solution of the master equation given by: 



Pit) 



/O 0\ 
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\0 0/ 



(16) 



This state is also an invariant state and its concur- 
rence is independent of time: C{p) = 1. 



Cl{p) = 2(^p23P32 - VP11P44) (13) 

C2{p) = 2(v'y^Ii^ - VP22P33)- (14) 



III. SOLUTIONS FOR INITIAL STATES IN DFS 

In this and the next two sections, all density matrices 
and expressions of concurrences will be referred to the 

{|0i), 102), 103), 104)} basis. 

a) Consider |0i) as the initial state. The solution of 
master equation is given by: 



Pit) 



/I o^ 
0000 
0000 

\o oy 



(15) 



This corresponds to an invariant state, and its con- 
currence is given by: C{p) = '^^2N^i ~^ which is a 
constant in time. 

The concurrence only depends of A^. For N = 
we have a factorized state at all times, but as we 
increase N , we get a maximally entangled state in 
the large N limit. (see, fig.l). 



In the following sections we consider, as initial states, 
103) OT 104) , and also superpositions of the form e 
|0i) + Vl - s'^ 104) and s |02) + Vl - s'^ |03). The idea 
is to increase s and to study the effect of having an in- 
creased component in the DFS on the death time of the 
entanglement. For simplicity, we assume 7 = 1, and 
i; = 0. 



IV. SOLUTIONS FOR A = 

a) The third initial state considered is |03). Initially 
its concurrence is: C(03(O)) = 1. It corresponds to 
a maximally entangled state. 

The solution of master equation for this initial con- 
dition and N = is given by: 



pit) = 



V 



,2t 



-l)e-2t 



e-2t 





(17) 



Since the matrix p{t)p{t) has only one nonzero 
eigenvalue, in this case we use the separability crite- 
rion [14]. According to this criterion, the necessary 
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condition for separability is that a matrix p^"^, ob- 
tained by partial transposition of p, should have 
only non- negative eigenvalues. In this particular 
case, we observe a negative eigenvalue for all times, 
so the state stays entangled, (Fig. 2). 
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which is shown in fig. 3 for various values of s: 
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FIG. 3: Time evolution of the concurrence for initial 
with : e = 0.28 (solid line), e = 0.345 (dashed line), e = 0.9 
(dotted line). 



FIG. 2: Negative eigenvalue (A) of p^^, for the separability 
criterion, for l^s) as the initial state. This eigenvalue is always 
negative , indicating entanglement at all times. 



b) Consider the initial state 104). Since = 0, 104) = 
I -h+) and C(04(O)) = 0, since is a factorized state. 

The solution of master equation for this initial con- 
dition is given by: 

^(-1 -2t-he2*)e-2* 



2te-2* 





Pit) 







-2t 



(18) 



and its concurrence is C{p{t)) = 0. 

c) Now, we consider an initial superposition of |0i) 
and 104) : |^i) = e |0i) + |04). 

As we increase starting from e = 0, we increase 
the initial projection onto the DFS. For s = 1 the 
initial state is in the DFS plane. 

For TV = we have: |^i) = s \ )-hVl -^^ |++). 

Its initial concurrence is C(^i(0)) = 26:^1 — s'^. 

The solution of master equation for this initial con- 
dition with = is given by: 



/ (-2t-l+2t£^+£^+e'^*) 

p2t 



Pit) 



\ 












Q 2t(l-£") 





e2t 









p2t 



J 



(19) 

and the corresponding concurrence is given by: 
C(p) =max{0, 2^14-/933} (20) 
= max{0, 2{{eVl-e^)e-' - te-^\l - e^))} (21) 



For 5 = and e = 1 , its concurrence is zero, there- 
fore we have a non-entangled state. For e > the 
initial entanglement decreases in time, and the sys- 
tem becomes disentangled (sudden death) at a time 
satisfying the relation: 



(22) 



For < £ < 0.34525 the equation (22) has two so- 
lutions, namely, when the system becomes sepa- 
rable, and tr > td when the entanglement revives. 
It should be noted that there is a critical e for which 
td — tr- For 0.34525 < e < \ the above equa- 
tion has no solution and the concurrence vanishes 
asymptotically in time. 

So, when we are "not far" from |04) we observe a 
sudden death and revival, but when we get "near" 
|0i) this phenomenon disappears. Fig. 4 shows the 
behavior of the death and revival time as function 
of e. 




FIG. 4: a) The death time b)The revival time of the entan- 
glement as a function of with initial |^i). 
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d) Finally, we consider an initial superposition of |02) 
and 103) : |^2(0)) = e |(/>2) + - l^s), which is 
independent of TV. So, like in the pervious cases, 
as we increase e, starting from 6 = 0, we increase 
the initial projection onto the DFS. For e = 1 the 
initial state is in the DFS plane. 

For = we have: |^2) = ^^[(^ + e'^) 

I — +) — (e — ^/Y^^^) I + — )] . Its initial concurrence 
is C(^2(0)) = |262-l|. 

The solution of master equation for this initial con- 
dition is given by: 



pit) 



\ 



e2t 





























(23) 



and the corresponding concurrence is: 

C{p) = max{0, |/)33 - P22I} 

= max{0, e-'^^\e^e^^ - 1 + 
which is shown in fig. 5. 



(24) 
(25) 



C(t) 




FIG. 5: Time evolution of the concurrence with intial \^2) 
, for : e = 0.3 (soUd Una), e = 0.5 (dotted hne), e = 0.707 
(dashed hne), e = 0.9 (dash dotted line) 



For < e < 0.707 the initial entanglement de- 
creases in time, and the system becomes disentan- 
gled instantaneously, at a time given by: 



1, 



(26) 



However, at the same time, the entanglement 
revives reaching asymptotically its stationary 
value, (fig. 6) 

When we approach the decoherence free subspace 
this phenomenon disappears. 

Next, we treat the cases with N > 0. 
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FIG. 6: Death-revival time as given by eq 26 , versus £. 



V. SOLUTIONS FOR iV ^ 



In general the solution of master equation for , 
104) and l^i) as initial condition and N ^ 0, has the 
following form: 



Pit) 



fpii pi4> 



P33 

Vp41 P44y 



and written at the standard basis, has the same struc- 
ture as in (12). Its concurrence is given by: C{p) = 
max{0, Cl(p), C2(p)}, with the explicit expressions for 
Cl{p) and C2{p) are given by: 



Clip) = 2( 



P33 



1 ipuN + p44iN + 1) 


+ 2pU^N{N + 1)) 


2N 


+ 1 


1 iPiiN + pii{N + 1) 


- 2pU^/N{N + 1)) 



C2{p) = 2( 



2N+1 

(27) 

WN{N + l){pn-Pi4) + Pu\ P33. 

2nT~i 



On the other hand, for the initial condition |^2) , the 
corresponding expressions for the density matrix and con- 
currence are: 



Pit) = 



/pii pi4\ 

P22 P23 

P32 P33 

\P41 P4J 
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CHp) = |P33 -P22I 

_ / N{pn + P44) + P44 + 2puyjN{N + 1) 
V 2N+1 

lN{pii+p44) + Pii-2pi4^/N{N + l) 
V 27V + 1 

(29) 

= ^^^|^iv(iv + i)(m - m) + P14I 

- a/ (/O22 - 2/523 + PZZ){P22 + 2/923 + Pas) (30) 

a) Next, we consider again the case for initial 
103) , but for N 7^ 0. Tlie concurrence is: 
C{p) = ■max{Q, C\{p), C2(p)}, with its initial value 
C(<^3(0)) = 1. 

In fig. 7, we show C{t) versus time for various val- 
ues of N . We observe sudden death in a finite 
time, then the concurrence remains zero for a pe- 
riod of time until the entanglement revives, and 
the concurrence reaches asymptotically its station- 
ary value. Notice that this time period increases 
with N. 

1,0-1 1 




t 

FIG. 7: Time evolution of concurrence for initial l^s) , with: 
N=0.1 (dashed line), N=0.5 (solid line), N=l (dotted line). 



In the fig. 8 we show the death and revival times 
versus N. They decrease and increase with N re- 
spectively. 

b) Consider 1^4) as an initial state. The concur- 
rence, in this case, is: C{p) = max{0, (72 (p)}, with 
C2{p) defined by (28). Initially, it takes the value 

C(04(O)) = '^^^j^^i'^ -- The behavior of concur- 
rence is similar as in 1^3). The initial entanglement 
quickly decays to zero, getting disentanglement for 




N 



FIG. 8: a) Death time b) Revival time versus N for the initial 
state I 03). 



a finite time interval , then the entanglement re- 
vives and asymptotically it reaches its stationary 
value. 

However, unlike the case with initial state 103), 
the death time first increases reaching a critical 
TV = 0.421, then decreases, as shown in Fig 9. The 
revival time has the same behavior as in 103). 




N N 



FIG. 9: Death (a) and Revival (b) times versus N for the 
initial state | 04). 



c) In the following case, we consider the superposi- 
tion l^i(t)) = e\(j)i) + \/l — ^^104) as the initial 
state. The solution of master equation for this ini- 
tial condition depends on e and N and also its con- 
currence, which is: C{p) = max{0, C2(t)}, where 
C2{t) is in 28. Since Cl{t) is always negative, 
the only contribution to the concurrence comes 
from C2(t), . Its initial value being C(^i(0)) = 
\2eVi £ +^V^(^+ i)^f — 2)1 ^ hence is it clear that 
for certain pairs of N and 6, our initial state will 
be a non-separable one. In the Fig. 10 we show the 
time evolution of the concurrence for N = 0.1 and 
several values of e. For e = and e = 1 we retrieve 
104) and |0i) respectively. For < 6 < 0.5 the con- 
currence dies in a finite time, stays zero for a time 
interval and subsequently revives, going asymptot- 
ically to its stationary value. For values larger than 
5 = 0.5, there is no more sudden death, since we are 
getting "close" to the DFS, and C{t) goes asymp- 
totically to its stationary value. The Figure 11 
shows the death times versus s foi N = {0, 0.1, 0.2} 
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FIG. 10: Time evolution of Concurrence for as initial 

state and A^ = 0.1: 5 = 0.1 (long dashed line), s = 0.2 (dash 
dotted line), e — 0.29 (dashed line), £ = 0.5 (dotted line), 
e = 0.9 (solid line). 



. There is a curious effect, that for Ny^O, as we in- 
crease £, the death time first decreases and subse- 
quently it behaves "normally", by increasing with 
e. In the Fig. 12 we show the revival time as a 
function of e for the same values oi N . In all cases 
the revival time decreases with e. 
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FIG. 11: Death time with initial and :Ar = (solid line), 
N = 0.1 (dotted Une), N = 0.2 (dashed line) 



d) Finally, we consider the case with initial |^2(^)) = 
e\(\)2) + Vl — ^^1^3) • Its concurrence is: C{p{t)) = 
max{0, Cl(t), C2(t)}, with Cl{t) and C2{t) defined 
in (29,30), audits initial value:C(^2(0)) = \2e^-l\. 
Fig. 13. shows the time evolution of the concur- 
rence with = 0.1 for several values of e . 

As we can see from the Fig. 13, this case is more 
complex, since there are more than one death and 
revival before reaching the critical value of e. Such 
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8 

FIG. 12: Revival time with initial |^i) and :iV = (solid 
line), N = 0.1 (dotted line), N = 0.2 (dashed line) 
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FIG. 13: Time evolution of concurrence for initial 1^2) and: 
£ = 0.1 (solid line), e = 0.4 (dotted line), e = 0.49 (dashed 
line), e — 0.54 (dash dotted line), e — 0.6 (long dashed line), 
£ = 0.9 ( space dashed line). 



a situation has been described previously [15, 17]. 
Like in the previous cases above a certain critical £, 
when we get "close" to the DFS, these effects disap- 
pear and C{t) goes asymptotically to its stationary 
value. 



VI. DISCUSSION 

The first and most obvious observation is that if we 
start with an initial state that is in the DFS plane, the 
local and non-local coherences are not affected by the 
environment, thus it experiences no decoherence and the 
concurrence is constant in time. It does increase with the 
squeeze parameter N (in the case of initial |(/>i)), getting 
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maximum entanglement for N-^ oo. So this reservoir is 
not acting as a thermal one, in the sense that introduces 
randomness. On the contrary, a common squeezed bath 
tends to enhance the entanglement, as we increase the 
parameter N. 

This is clear if we observe that for oo, 

^(| + +) + | )), which is a Bell state. On the other 

hand, if we start with the initial state 102), this state is 
independent of N and it is also maximally entangled, so 
C=l for all times and all N's. 

Now, we consider other situations with initial states 
outside the DFS: 

For N=0, |0i) ^ I - -); 104) ^ I + +) while |02) and 
103) are independent of N and represent the antisymmet- 
rical singlet and one of the triplet states, respectively. 

If our initial state is |03), the initial concurrence is 
C(t=0)=l, and for the steady state, C(t=oc)=0. On 
the other hand, from the negativity of the eigenvalues of 
p^^ (Peres-Horodecki criteria), we know that the state is 
entangled for all times, so presumably the concurrence 
versus time is a smooth curve starting from 1 and going 
asymptotically to zero. 

For an initial |04) = | + +), of course C(t=0)=0, since 
we have a factorized state. Since this initial state is sym- 
metrical, the steady state is |0i) = | ), which is also 

a factorized state with C(t=oo)=0. Our analysis shows 
that C(t)=0 for all times. 

Finally, if we have initial states of the form |^i) = 
5|0i) + ^1^104 ) = e I - -) + I + + ) or 

1^2) = s|02) + VT^|03) = ^[(e + Vl^) I - 

-\-) — [e — \/l — I H )], there is sudden death if 

< 5 < 0.34525 and < 5 < 0.707 respectively. Obvi- 
ously, when we get "near" the DFS, that is e gets larger 
than zero, the sudden death times become larger up to 
some critical value, for which the death time becomes in- 
finite and we no longer observe the sudden death. Thus, 
we have related directly the local decoherence with the 
disentanglement . 

Another interesting feature is the revival. For the |^i) 
case, after a finite time period, the entanglement "re- 
vives", and the concurrence reaches asymptotically its 
steady state value. This time gap becomes smaller for 
larger e until we reach e = 0.34525 . At this point, the 
two times merge and beyond that, there is no longer sud- 
den death nor revival. 

On the other hand, for the 1^2) case, the sudden death 
and revival happen simultaneously, thus the above men- 
tioned time period vanishes. The phenomena of one or 
periodical revivals have been obtained before, but always 
in the context of one single reservoir connecting both 
atoms, like in the present case [15-17]. 

In the Nt^ case, we also observe sudden death and 
revivals, up to a certain "distance" to the DFS ( or more 
precisely, up to a certain critical value of e). 

However, there are certain differences with the N=0 
case. 



For the initial state |04), the death time versus N first 
increases for small values of N, and for N> 0.421, it tends 
to decrease. A possible interpretation of the increase is 
the following one: 

l^4) = 7^i^(M| + +) - 7V| - -)) = 

Vivw (7^ + ^l + +)-!- -)), so that the ra- 
tio of the probabilities of the double excited and the 
ground states goes as (1 + -^.) 

On the other hand, the squeezed vacuum has only com- 
ponents for the even number of photons, so the interac- 
tion between our system and the reservoir goes by pairs 
of photons. Now, for very small N, the average pho- 
ton number is also small, so the predominant interaction 
with the reservoir will be the doubly excited state that 
would tend to decay via two photon spontaneous emis- 
sion.. Now in the |04) case, the population of the | + +) 
goes down with N, meaning that the interaction with 
the reservoir goes also down with N and therefore, the 
death time will necessarily increase with N, which de- 
scribes qualitatively the first part of the curve (fig 9-a). 
On the other hand, as we increase the average photon 
number N, other processes like the two photon absorp- 
tion will be favored, and since there will be more photons 

and the | ) population tends to increase with N, this 

will enhance the system-bath interaction and therefore 
the death of the entanglement will occur faster, or the 
death time will decrease. 

In the 103) case, initially there is no | + +) component, 
thus we expect a higher initial death time. However this 
case is different from the previous one in the sense that 
the state is independent of N, so there is no initial in- 
crease. However, as the state evolves in time, | + +) and 

I ) components will build up and the argument for 

the decrease of the death time with N follows the same 
logic as in the previous case, (fig 8-a). 
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VII. APPENDIX 

In the representation spanned by 
{101)7 102)5 103)5 104)}, the solution of master equa- 
tion (6) is: 
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P23W=P23(0)e-(2^+l)* 



Pll{t) = 



^N{N+1){8N + 4) 



{[-(2Ar+l)(p44(0)+p33(0)) 



P24W=P24(0)e-(2^+l)* 



(P44(0) + P33(0))Ar + ip44(0)]e-^*(^+^)' 



-8e-^*e-(2JV+i)t 



+ [-(2Ar+l)(p44(0) + p33(0)) 



A/'{A/' + l) 



N{24:N3 + 36Af2 _,_ iqa^ - 1) 
^V2N + W2N^ + N)(N + ^) 



[{nVn + T 



{P44{0) + P33m)N^ - (/944(0) + P33{0))N 



-P44(0)]e 



+ {2N + l)(p44(0) + P33(0) + pn{0))VN{N+l)} 



Pl2(t)=Pl2(0) 



12e-(2^+i)' 2 1 

iV(24iV3 + 36iV2 + lOA/" - 1) 3^ ^2' 



-(((2A^+l)t + J1')V2jV2 + iV + 4t7Vv/iV+lV2iV + l 

X (e^*p34(0) + P43(0))e 

+ e**(-;V + ]^){N^/WTl + ^ V2Ar + l\/2N'^ + AT) 

-(((2Ar+l)t + J^)\/2]v2 + jV-4tArViV + l V2iV + l 

X (e^*/534(0) - P43(0))e V^^^ 

- ^ vW+lV2iV2+iV(p3i (0) e'* (TV^ + iV - i ) 

- \{N + \)pim) - m + l)iVV]vTlp34(0)e^*] 



X {NVnTT+ ^y'2N + l\/2N^ + N) 

^ (-^y2N'^ + N(2N+l)+4Ny/N+ly/2N+l) 

x(/,43(0)-e**/>34(0))e 

- |(iV + ^){NVn + 1 - ^V2N + W2N^ + N) 

^ {^y2N'^ + N{2N+l)+4:N^/N+l^/2N+l) 
y/2N'^+N 



P32W=P32(0)e-(2A^+i)* 



/'33 



X (e'*p34(0)+p43(0))e" 



(*) = ^[e-^*(^-^)^(p33(0)+P44(0) 
+ [e-2*(^+^)'(p33(0)-p44(0) 



Pii{i) 



+ ^)p34(0)e^* + pi3(0)(Ar2 + TV - ^)) 
^/2N + lv/2Ar2 + TV + ^(AT + i)V]vTl/943(0)] 



gg-(2iV+l)t 



v/Af(Ar + i)^ 



= ^(^34(0) - e-^*/,43(0))e-2*(^-^)^ 
+ ^(P34(0) + e-*p43(0))e-2*(^^+v^)^ 



(2Ar + i)(i2Ar2 + i2Ar- 1) 



jg-(2JV+l+4y'jV(JV+l))i 



X {-i^{N + 2)(2/O44(0) + P33(0))\/iV(iVTT) 

+ (i/944(0) + P33(0))(A^2 ^ ^) ^ l^^^^o^^ 
_ g-(2JV+l-4VAr(JV+l))t(-l(^^ _^ 1 ^ 

X (2/>44(0) + /933(O))0V(A^fl) 

+ (^P44(0) + P33(0))(iV2 ^ ^) ^ \pim) 

+ ^(V2]vTlpi4{0){iV2+iV- ^)V2A^2 + jv 



P4l(t) = 



12e-(2^+i)* 



A^\/]VTl(2A^ + l)(12Ar2 _^ i2A^ - 1) '3 



[-7V^/]vTl 



+ -{N + -)(2p44(0) + P33(0))ArA/iV + l)] 



P2l{t)=P2m 



P22{t) = P22{0) 



X (-(27V + 1)(^ + p^mVNiN + 1) 

+ (2p33(0) + /944(0))(iV2+iV) 

P44(0) X (2Ar+i+4^/jV(jV+l))t _ Ij^^J^ , 1 

4 ^ 3 
X ((2iV + 1)(^ + P4i{0))VN{N + l) 

+ {2p33{0) + p^miN^ + N) 

_^ P44(0) .j^^-(2N+l-4^/N(N+l))t 

4 

+ v/iV(ivTi)(\/2ivTip4i(o)(iv2 - ^ + ^) 

X V2iV2 + N+ -{N + -)NVN^N 

X(^+P44(0)))] 
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P42W=P42(0)e-(^^+l)^ 

+ ^(P43(0) + e^^p34(0))e-^^(v^+^)^) 
WO) V7V(iV + l) 
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